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Abstract: The most general spherically symmetric solution with zero shift is found in the non- 
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conventional thermodynamics of this general class of solutions is established by calculating the energy, 
temperature and entropy. Although several of the solutions have conventional horizons, for particles 
with ultra-luminal dispersion relations such solutions appear to be horizonless. 
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1. Introduction and results 

The UV completion of gravity has been a difficult road for theoretical physics of the past fifty 
years. The only convincing answer so far has been provided by string theory, and it works only in 
perturbation theory and at energies well below the Planck scale. 

Recently, a different field theory model for a UV complete theory of gravity was proposed ||1[] 
(see also 0). The theory does not have the full diffeomorphism invariance of GR but only a subset. 
Because of this property a different scaling symmetry is allowed in the UV and the theory accepts 
renormalizable couplings with up to six derivatives. It is interesting that this picture allows a theory 
of gravitation that is scale-invariant in the UV, and where standard general relativity with its higher 
symmetry could be an emerging theory in the IR. 

There are several versions by now of the Hofava-Lifshitz type of theory. Two main categories 
were introduced originally by Hofava. The non-projectable version allows the lapse to be a general 
function of spacetime coordinates. In the projectable version, is a function of time only. 

Furthermore in the original formulation the principle of detailed balance (DB) has been imposed. 
It had the advantage of reducing the possible coupling constants of the theory. In subsequent works 
[D, ^, it was advocated that a general action should be used that allows all couplings compatible with 
renormalizability. It was ffist argued in that DB implies a correlation between the effective Planck 
scale and cosmological constant that left the theory very little room to agree with observation. 

It was subsequently shown in 0, that the large distance gravitational field of a spherically 
symmetric source has a very different behavior in the HL theory with detailed balance than GR, a 
fact that excludes HL gravity with detailed balance as a description of low energy observable gravity. 
At the same time, a small breaking of the detailed balance reinstated the usual 1/r tail of gravity ^j. 

It is obvious that we have a family of HL-type theories with projectable or non-projectable lapse, 
and with or without detailed balance. In view of this, we will reiterate the philosophy advocated in 
0: the most general HL theory compatible with renormalizability should be investigated, not only 
in terms of its consistency but also in terms of its ability to agree with observables in gravity and 
cosmology. 

The arena where the theory seems most promising is cosmology. It was initially pointed out that 
the UV structure of the theory can help in solving the horizon and flatness problems without inflation 
0. In particular, the fact that speed of light in the UV is infinite indicates there is no horizon problem. 
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while the fact that curvature contributions to the Friedmann equation are enhanced, indicates that 
they must be suppressed at early times. This is similar to standard holographic cosmology 0, 
although the mechanism for the generation of the curvature square corrections is different. 

Scale invariant perturbations can be generated naturally in the HL high-energy phase without 
inflation 0,0. Relevant terms provide power suppressed corrections to the scale invariant spectrum. 
It was argued in that the logarithmic running of marginal couplings in the UV will provide a 
logarithmic tilt in the scale invariant spectrum of perturbations, although the direction of the tilt 
cannot be determined without a detailed calculation. 

The presence of higher derivative terms in the theory provide the opportunity of a bouncing 
universe 0, 0], [|^. The existence of parity odd couplings at high energy provide the possibility of 



polarization asymmetry in the CMB data |]rT|. It should be noted that CP non-invariance is not a 
necessary condition of the formalism, but it is certainly allowed. It is present in the original Hofava 
proposal. 

In the projectable theory, the Hamiltonian constraint is not a local constraint. It was shown in 
[0 that in the cosmological context the constraint allows dark matter as an integration constant in 
the FRW equations. Many further analyses of cosmological issues have been done since ||13|| . 

The issue of spherically symmetric solutions and black holes, as well as issues related to thermo- 



dynamics, were analyzed in several works |14H16]. Other issues involved strong coupling problems 



[0-1231, comparison with solar system data ||2^, issues of quantum field theory and renormalizability 
[ p^ , as well as the discussion of particle geodesies [p^-p8|. 

An important issue is related to the reduced diffeomorphism invariance and the associated strong 
coupling problems. As anticipated in [Q, the dynamics of the extra scalar mode that exists in the 
theory is particular and can lead to strong coupling problems. This was argued to be the case first 
in [TH] and most convincingly in [BO]. Moreover, generic perturbative instabilities were found in the 



renormalizable regime of couplings p2|. The projectable theory seems to be free of such problems, 
but is expected to have problems with caustics in the cosmological domain |2^. In arguments 
against caustics formation were advanced. Recently a modified theory, where additions of spatial 
derivatives of the lapse were utilized, was argued to be problem-free . 

An important issue for this class of gravity theories remains open. This is the issue of quantum 
UV structure and renormalizability. Although the theory is power-counting renormalizable, several 
of the important UV couplings are dimensionless. They will therefore have logarithmic UV running, 
and renormalizability requires that such couplings are asymptotically free. Although this can always 
be guaranteed by taking the appropriate sign for a given coupling, positivity of such couplings may 
not allow the appropriate choice. For example, A > 1 for the phenomenological viability of the 
theory. The functions for the UV marginal couplings have not yet been calculated and therefore 
this important issue remains unresolved. In particular, the question whether A = 1 is an IR fixed 
point remains open. 

The potential asymptotic freedom of the UV couplings of the theory may have an additional 
impact on the cosmological constant problem. Indeed, it allows the possibility that an exponentially 
small vacuum energy is generated because of asymptotic freedom. 

An important coupling in this theory is the speed of light. This is a relevant coupling in the UV 
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theory and therefore its RG running will be powerlike. It is however a marginal coupling in the IR 
theory and therefore its IR asymptotics will have a logarithmic running. This implies that the RG 
behavior cannot be followed from the UV to the IR by perturbation theory. 

The breaking of Lorentz invariance in the gravitational theory leads to several important and in 
principle observable effects. As was pointed-out first in and subsequently investigated in many 
papers, the matter sector of the gravitational theory will have Lorentz invariance broken, either by 
the UV couplings or by quantum effects that communicate Lorentz invariance breaking from the 
gravitational sector. There are severe constraints on such Lorentz violating couplings that have been 
mostly derived in the last decade or so ||29|-^^. A particularly important question is to what extend 
the speeds of light relevant for different particles will be equal. Inequality of such speeds of light 
in the IR is strongly constrained. On the other hand, the RG structure of non-relativistic matter 



theories is very interesting |33], and may also accommodate novel effects for the SM |34 



Issues due to the broken Lorentz invariance arise also in the study of spherically symmetric solu- 
tions. Indeed, in such cases the notion of the horizon may be particle dependent if Lorentz invariance 
is broken. This is an issue that can be studied using the appropriate generalization of geodesic 
equations as in [^-[^. This however implies that the notion of temperature is therefore particle- 
dependent, if particles have dispersive geodesies. Worse, the notion of thermodynamics becomes fuzzy 
as the analogue of Hawking radiation depends on dispersion and ceases to be exactly thermal ||35| , 
3^ . Moreover instabilities appear through the form of perpetuum mobile |^ . 



It has been speculated by many, that non-relativistic gravitation theories of the Hofava-Lifshitz 
type may serve as duals of non-relativistic strongly coupled large -A^ quantum field theories. We 
find such expectations to be remote and only for field theories that break translation invariance. 
Standard holography as described by standard (super-gravity) has all the ingredients to describe 
translationally invariant non-relativistic quantum field theories of the Lifshitz type. The reason is 
that such QFTs are different from Lorentz invariant ones by the choice of couplings in the Lagrangian. 
Indeed, by choosing non-Lorentz invariant sources breaking the symmetry between space and time, 
we can generate in the UV the appropriate non-relativistic QFT. In the dual (supergravity) language 
the equations remain the same (fully diffeomorphism invariant), only the sources change. 

Another hint in the same direction is provided by the intuition that it is the translation invariance 
of a boundary QFT that is promoted to diffeomorphism invariance in the bulk leading to standard 
gravitational bulk theories. As the diffeomorphism invariance of Hofava-Lifshitz gravities is smaller, 
this suggests that such a gravitational theory maybe dual to a QFT with reduced translation in- 
variance. A better understanding of the holographic role of Hofava-Lifshitz gravity remains to be 
found. 

The purpose of this paper is to find and analyze spherically symmetric, static solutions in Hofava- 
Lifshitz gravity. The appropriate ansatz for the metric degrees of freedom can be written as 



where fi^ is a two sphere of unit radius when = 1, a torus of unit volume when k = and a 
pseudosphere of unit radius when k = —1. Unlike the fully diffeomorphic invariant case, we cannot 



-4 - 



set the shift A^^' to zero by a coordinate transformation. 

The general action of the Hofava-Lifshitz class of theories is^ 



S = j dtd^x^N[a{KijK'^-XK^)+^£'^''RuVjR\+CRijR'^+vR^+^R+(^]+S3, (1.2) 

where S3 is the part of the action with six derivatives 

83 = j dt£x^N [pCijC'^ + piRDR + p2R^ + p3RRijR"^ + p4RijR"''Rh] ■ (1-3) 

We will make the simplifying assumption in this paper that f3i = (^2 = = = 0. The reason is 
that apart from simplifying the equations, the qualitative structure of the solutions we will find is 
expected to be similar. Indeed, the presence of all curvature square terms takes into account curvature 
non-linearities. 

To make the problem tractable we will also look for solutions that have A^^ = 0. In this case we 
will find the most general solutions to the non-linear equations of the theory. 
We find that: 

• It is a generic feature of the solutions found that they have regular large distance asymptotics 
that are asymptotically AdS/dS or flat. Moreover, generically the next correction is compatible with 
a standard Newton's law. 

There are exceptions to this result (in special cases). In particular, the detailed balance action 
first written down by Horava is one of these notable exceptions. As shown already first in 0, it 
does not reproduce the correct Newton law at large distances. This is part of a class of special cases 
analyzed here, occurring when 

B = ^(^-2.)=0 or C = 15^^<0. (1.4) 

although in the second case it is not the generic solution. 

• For the general solution we calculate the temperature, mass and entropy, assuming Lorentz- 
invariant probe particles with standard dispersion relations and using the first law. In one of the 
two general categories of parameters {A < with A defined in ( |6.2| )) where the horizon distance 
is permitted to shrink to zero, the entropy has a logarithmic divergence in this limit. In the other 
general category {A > 0) the horizon distance is bounded from below. 

• For special values of the curvature- squared couplings ( + 3rj = 0, when the horizon position is 
allowed to go to zero, the entropy is regular in this limit. 

• For the same values of the curvature-squared couplings C+3?7 = 0, there exists a solution with 
logarithmically corrected large distance asymptotics. Its "mass parameter" defined in the naive way, 
depends logarithmically on radial distance and becomes smaller at larger distances. 



^This general action was advocated in on the basis of analyzing the phenomenology of Hofava-Lifshitz gravity. 
In /3i = as contributions to FRW equations were studied. /3i ^ was introduced in [Q where a different basis of 
marginal terms is used. 
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• When the cosmological constant term in the action is absent, the effective cosmo logical constant 
of (A)dS curvature is also zero (for some branches). 

• We study geodesies of particles with finite and infinite light speeds in spherically symmetric 
backgrounds with traditional horizons. We parameterize the dispersion relations of particles as pg = 
(p^)"", with n > 1. We find that when n > 1 the traditional behavior of the horizon disappears, 
suggesting that for such particles the black hole is effectively "naked" . 

The structure of this paper is as follows. 

In section 2 we review the Hofava-Lifshitz theory and its generalizations. 

In section 3 we introduce the ansatz for the solutions, as well as the equations to be solved. We 
also include a discussion of several important conceptual issues for such solutions ranging from the 
notion of horizon and singularity to the notion of black-hole thermodynamics. 

In section 4 we discuss black hole solutions in the standard case, where higher curvature correc- 
tions are set to zero in order to establish notation and connect to the standard cases. 

In section 5 we discuss the spherically symmetric solutions in two special cases, determined by 
specific relations of the curvature-squared couplings. In the first case our solutions are new. In the 
second, our solutions generalize previous solutions found by Lu, Mei and Pope ||^, Kehagias and 
Sfetsos 0, and Park [|l5l. 

In section 6 we discuss the most general spherically symmetric solutions, by distinguishing two 
main cases involving the coupling constants. 

In section 7 we describe the analysis of the conventional horizons for the general solutions. 

In section 8 we calculate the conventional thermodynamic properties (for regular probes) assuming 
the validity of the first law. 

Finally, section 9 contains a summary of the results as well as open questions. 

In four appendices we collect several technical calculations. 

Appendix A contains a derivation of the general equations of motion. 

Appendix B fills-in the details for the derivation of the special class of solutions. 

Appendix C fills-in the details for the derivation of the general class of solutions. 

Finally, appendix D derives "generalized geodesic equations" and solves for the radial geodesies. 



2. The Horava-Lifshitz gravity theory and generaUzations 

We review here the Hofava-Lifshitz gravitational theory as was formulated in and as it was 
generalized in subsequent works 

The dynamical variables are N,Ni,gij, with scaling dimension zero, except Ni that has scaling 
dimension 2. This is similar to the ADM decomposition of the metric in standard general relativity, 
where the metric is written as 

rfs2 ^ _j^^2 ^ ^ N'dt){dx^ + N^dt) , Ni = QijN^. (2.1) 

The scaling transformation of the coordinates is now modified to 

t^et , x' ^lx\ (2.2) 
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under which Qij and are invariant, while scales as A^* — > £ ^A^j. 
The kinetic terms are given by 



SK = ^j dtSx^N {K,,K'^ - \K^) , K = g'^K,, , K'^ = g'^g^'Ku (2.3) 



in terms of the extrinsic curvature 



K.J = ^{g.j - ViN, - VjNi), (2.4) 



with covariant derivatives defined with respect to the spatial metric gij. 
The most general power-counting renormalizable action is 

S = j dtd^Xy/gN[a{KijK'^ -XK^)+-fS'^''RiiVjR^k+CRijR"^+vR'^+^R+(^] + (2.5) 

with S^^'' = ^ the standard generally covariant antisymmetric tensor, e^^^ is defined to be 1, and 
other components are obtained by antisymmetry. Indices are raised and lowered with the metric gij. 
Therefore, S^^'' = (±l)/y^. 5*3 is the part of the action with six derivatives 

S3 = j dtd^x^N [PQjC'^ + piRDR + p2R^ + PsRRijR'^ + p4RijR'''Rh] ■ (2.6) 
The action with detailed balance corresponds to the following values for the coefficients above 

''^^=8(1-3A) 4 ' ^^^^S^T^"^^ ' '''''' = W^y '^^ ^ ^ ^ 

as well as Pi = P2 = P3 = Pa = 0. 

The action in (|2.11|) is invariant under a restricted class of diffeomorphisms 

t' = h{t) , x'' = h\t,x). (2.8) 

The transformation of the metric under infinitesimal diffeomorphisms is 

Sgij = die^gjk + dje^gik + e^dkgij + fgij (2.9) 

6N, = d.e^N^ + d.e^N, + g,, + fN, + /iV, , 6N = eW,N + fN + fN. (2. 10) 

This is the first version of the Hofava-Lifshitz theory we will analyze in this paper. A second 
version, called projectable, assumes that A^ is a function of time only and can therefore be set to 
one (if non-zero) by a diffeomorphism. We will discuss spherical solutions to this second theory in a 
subsequent publication. 

To proceed further with our solutions we will make the simplifying assumption that Pi = P2 = 
Pz = Pa = 0. The reason is that apart from simplifying the equations, the qualitative structure of the 
solutions we will find is expected to be similar. Therefore, the action we will consider is 

S= j dtd^Xy^N[a{KijK'^-\K^)+pCijC'^+-f£'^^RuVjR^k+CRijR'^+vR^+^R+(^], (2-11) 
where the various coefficients are assumed independent. 
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2.1 The IR limit around flat space 

Around flat space the IR action simplifies 

S ^Se = j dtd^x^N [a{KijK'^ - \K^) + ^i? + a] . (2.12) 
Defining = ct, choosing A = 1 and 

c=\ - , MI^IQitG = ^ , A = -;^, (2.13) 
the action is that of Einstein 

= I d'xy^ [K,^K'^ -K' + R-2A]= -i- / d^xy^ \r - 2a\ . (2.14) 

Io/tG J IottG J L J 

The full space-time metric (jfj^i, is 

goo = -N^ + Nig'^N, , go^ = N, , gij = Qij , det[^] = det[^]iV2^ (2.15) 
while the inverse metric g'^" 

Q = ^ , Q = — n , Q = Q (2.16) 

3. Static, spherically symmetric solutions with zero shift 

The most general static and spherically symmetric metric can be written as 

ds^ = -(N(rf - NJrf )fir)de + 2NJr)drdt + -77- + r^dQl, (3.1) 

fir) 

where dfll is the metric of a two-dimensional maximally symmetric space. For /c = 1 it is a sphere 
of radius 1, for = — 1 it is pseudo-sphere of radius one, and for /c = it is a two-dimensional torus. 
We have also parameterized the lapse A^^ = N'^f for further convenience. 

Unlike standard general relativity, the reduced diffeomorphism invariance here is not enough to 
set Nr to zero. However to simplify the equations we will consider the ansatz where A^^ = 0. This 
will give as only a subset of all possible solutions. The metrics therefore we consider are 

dv'^ 

ds^ = -N{rff{r)de + 7^ + r'dnl (3.2) 

We now substitute this ansatz into the equations of motion for the (source-less) action (p.ll|) . They 
were first derived in 0] and reproduced here in appendix ^ 

Note that since for the metrics (|3.2| ) both Kij and the Cotton tensor vanish, it is only the following 
part of the action that contributes to the non-trivial equations 

S=j dtd^Xy/gN [CRijR'^+r]R'^+^R+a] , (3.3) 
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Therefore, only the coefficients C^ri,^,a enter the relevant equations. 

We will assume that ^ > so that we have regular GR at weak curvatures. Because of the 
existence of the couplings cubic in the curvatures, there are no obvious constraints on the signs of 
C,r). Subtler constraints will be obtained later from the behavior of non-trivial/non-linear solutions. 
We will be mostly interested in solutions that have either asymptotically AdS or asymptotically 
flat behavior, although our results are valid also in other cases. There are cases for example where 
the asymptotic behavior is dS, but in such cases the natural solutions should be time dependent. 
The asymptotic behavior is not only governed by the "cosmological" constant a but also by the 
curvature-squared couplings (,1] via the known phenomenon of self-acceleration. 

In the sequel we will use units of time so that c = 1 which via ( p. 131) implies a = ^. 

The N equation implies 

(3C + 8r/)rV^ + 4r(/ - k) [{C+Ar])f' - ^r] - A^r^f + 4(C + 2r/)(/ - k)^ + 2ar^ = 0. (3.4) 

The 12, 13 and 23 equations are trivially satisfied. The 11 component of the spatial (ij) equations 
gives 

^iiiV' + 5iiiV = 0, (3.5) 

where 

An = 4rf [{3C + 8r])rf + 2{C+4r]){f - k) - 2ir^] (3.6) 
5n = (3C+8r7)r2(4/r+/'2)+4(/-A:)[(C+477)r/'-er2-2A:(3C+8r7)]+2r=^((Tr-2e/')-4(5C+14r7)(/-A;)2. 

(3.7) 

The 22 equation (the 33 is identical to 22 due to symmetry) is not independent but can be derived 
from ( p.4|) , ( ^l5|) by differentiation. 

Using (|3^) , the expression ( p.7|) assumes the simpler form 

5n = 4(3C + 8r/)/ [r' f" - 2{f - k)], (3.8) 

and equation ( ^TSf ) becomes 

in(lniV)' + 5n = 0, (3.9) 

where 

in = r[(3C + 8r/)r/' + 2(C + 477)(/ - k) - 2^r'] (3.10) 

5n = (3C + 8r/)[rY'-2(/-fc)]. (3.11) 

Therefore, we end with two equations ( p.4|) and ( p.9|) for the two unknown functions f{r), N{r). 

In order to solve equations ( p.4|) , (|3.9|) , we will start with appropriate special cases in the space 
of parameters (, t], and will then proceed to the general case. 

3.1 Large distance asymptotics 

In most of the solutions we will find, we have recognizable large distance asymptotics of the form 

/(^) = A; - - ^ + 0(r-4) , = f = k - - + Oir~^). (3.12) 

3 r 3 r 

In all cases with such asymptotics M and M are the same. There are however special cases where 
the large distance asymptotics above are not valid. This included the case of detailed balance. 
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3.2 On horizons and singularities 

In standard general relativity and its generalizations that keep diffeomorphism invariance intact, 
horizons and singularities play an important role and their presence is relatively easy to discern. 
Although there are various types of singularities from serious to mild ones (where all curvature 
invariants are regular) the tools for their detection involve geodesies and curvature. 

In theories without full diffeomorphism invariance, both of the previous concepts are harder to 
discern. The fact that particles can have non-standard dispersion relations, and therefore no uniform 
maximal speed, implies that the notion of a horizon may be different, and indeed, a concept that 
is probe dependent. Indeed, in appendix we analyze geodesies of particles with non-standard 
dispersion relations, which suggest that in ultra-luminal cases there is no horizon in the standard 
sense. 

A different issue is to what extend "horizons" as defined by the divergence of metric components 
are indeed regular. In diffeomorphism invariant theories, a change of coordinates suffices to make 
the case. In theories with reduced diffeomorphism invariance, this may not be always the case. 
Indeed, several well behaved coordinate systems involve coordinate transformations that are not 
symmetries. Moreover, even coordinate systems (like the tortoise coordinate) that can be reached by 
defacto symmetries may be not well behaved coordinate systems. This implies that for particles with 
standard dispersion relations the horizon may be a singularity. 

Finally, there is very little known about the relation of singularities with curvatures in such 
theories. Indeed, in HL theories in particular, both 3 and 4-dimensional curvatures can be computed, 
but as the previous argument suggests, there may be more sources for singularities. 

It is fair to say we know little in this direction, but the problem is of fundamental importance. 
We will start its study in this paper although our results are so far modest. 

3.3 Thermodynamics 

Black-hole solutions with horizons in general relativity have a well understood thermodynamic for- 
mulation that involves the first law and the associated Bekenstein-Hawking entropy. The mass of 
the black hole can be computed from the asymptotic behavior of the solution. The temperature 
can be computed via the Gibbons-Hawking technique, by rotating to Euclidean space and asking for 
regularity at the horizon. Finally, the entropy can be computed from the horizon area. 

In the case of general non-relativistic black holes, the notion of thermodynamics variables are 
not obviously well-defined. To start with, the notion of a horizon depends on the probe particles 
involved. For example, if the theory contains two types of particles with different relevant "speeds of 
light", then each particle sees a different horizon. The upshot is that entropy, if defined in terms of 
the horizon area, is particle-dependent. 

There are more complicated cases where particles have deformed dispersion relations and some- 
time no upper limit on the speed. This indeed happens for example in the UV regime of the HL 
gravity theories. In such a case the notion of the horizon seems not that well-defined. However, in such 
theories the notion of a geodesic is modified and propagation should be computed from scratch. Such 
issues have been discussed in several articles |^ , , and we leave them for future exploration. 



- 10 - 



Not only the entropy but the notion of temperature is also probe dependent in the absence of 
relativistic invariance. In ||36| the Hawking radiation spectrum was computed assuming particles have 
modified dispersion relations and found to deviate from the thermal distribution. 

Finally, there is always the question of validity of the first law of thermodynamics, but experience 
so far seems to indicate that it is generically valid provided all ADM data are included. 

All of the above have been linked to wither strange or unacceptable behavior in non-relativistic 
black holes. In particular, in simple cases, it was shown that energy can be pumped out indefinitely, 
casting doubts on particular realizations p7| . 

All such issues need to be explored in depth, but we will not do this in the present paper. This 
is why we call the spherically symmetric solutions we found "black holes" , as it is not clear when and 
for whom they are black. 

In the sequel, we will "define" our thermodynamics by assuming that all relevant probes have 
standard Lorentz-invariant dispersion relations, as this guarantees that the temperature calculated a 
la Gibbons- Hawking is relevant^. Moreover, it will identify horizons as the largest root of the function 
/(r). This will coincide with the largest root of the function N{r). 

We will also assume that the first law of thermodynamics is correct. This was shown for spherically 
symmetric cases in ||3^ but the general case is still open. It is important to note that the proper 
definition of the thermodynamic energy is given by the mass parameter M in ( p.l2|) . The reason 
is that it scales properly under constant rescaling of the time coordinate. Knowing the mass and 
temperature allows us to compute the entropy up to an additive constant that we cannot fix from 
first principles. 

Note that in higher derivative theories of gravity that are fully diffeomorphism invariant we have 
the Wald formula for the entropy . Unfortunately, it is not at all obvious if such a formula or its 
modification is applicable in our case. 



4. No quadratic curvature terms 

• C = ^7 = 0. This is the simplest case, and reduces to Einstein gravity plus cosmological constant. 
Indeed, equation (|3^ ) becomes 

rf' + {f~k)-^r' = 0. (4.1) 

Its solution is 

^ , ^ 2 c , Ar2 2GM 

fir) = k + -y-- = k-- , 4.2 

0^ r 6 r 



with c an integration constant, the mass M = ^ , and G, A have been defined in ( |2.13| ). In particular, 
since we use units in which c = 1, we have that 

= IGttG = ^. (4.3) 



^This is obviously false with gravitons, and such thermodynamics is inapplicable in their case. 
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From (|3.11|) , we obtain Bn = 0, therefore, the solution of equation ( |3.9|) is N{r) 
Rescahng time we can set this constant to one so that 

N{r) = 1. 

This is the standard (A)dS-Schwarzschild solution. 

5. Special quadratic curvature combinations 

There are two possibilities that must be analyzed separately: 

5.1 C+377 = 0, 



No a constant. 



(4.4) 



2(T 

V 



In this case, equation ( |3.4| ) becomes 

.1] ■ J Li] 

We define a new function y and a new coordinate R 

2 A 
— r + 

6C ^ 
Equation ( |5.1[ ) becomes 

2/2 + 4y + 12y = 0, 

where a dot is a derivative with respect to R. 

The general solution of equation ( [5.3| ) is given by the implicit expression 

^ T \ 3 
' c 

r 



rV" + 4 V-(/-A;)lr/' + 4[^r2 + (/-A;) (/-/t) - — = 0. 



/(r) = A; H — -r^ + —r'^y(r) , i? = Inr . 



(5.1) 

(5.2) 
(5.3) 

(5.4) 
(5.5) 



where e = ±1 and is an integration constant (positive or negative). The modified lapse is 

ge(v^T=%-l) 

N{r) = , 

where we have fixed the multiplicative integration constant appropriately. The details of the derivation 
are given in appendix p. 

We will now study the large distance, r ^ oo behavior of the solutions above. 

• For large distances r » |ro|, the e = 1 branch of ( p. 41 ) can be approximated as y ~ r~^, and 
therefore 

/(^) ^ A; + — r2-^^ + 0(r-3) , GM = -^rl (5.6) 
6^ r 36 7] 

For an asymptotically AdS or flat solution we must therefore impose a > 0. 
For the modified lapse we obtain 



N{r) 



2e2r6 



O (r-9) . 



(5.7) 



This implies that here M = M. 
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y{r) ~ -3 
from which we obtain 



1, 


the laig 


:or 


To > 0) 






lo 









+ , N{r) 



log^ (r^/r^) I V r 



+ 



(5i 



f(r) ~ /c H H 



+ 



(5.9) 



Expanding the lapse at large distance we obtain 



A^^ ~ H + 



^ 2 , 



6^ Sr] 



log 



- log log 



2GM{r) 



+ C(r-^) (5.10) 



with 

2GM{r) 



elog^(r3/r^) 



cr 
6^ 



log^(r'^/r3) 



.(5.11) 



In this case there are logarithmic corrections to the effective cosmological constant as well as the 
Newtonian mass. They may be interpreted as the mass running logarithmically, and becoming 
smaller at large distances. We will not consider this solution further as its interpretation is not 
clear. 

In the above analysis we have assumed the coupling of the Einstein term ^ 7^ 0. If we evaluate 

(5.12) 



the solution for couplings relevant for detailed balance as in ( |2.7| ) with A = 00, we obtain 

CdB + ^DB = , ^DB = O-DB = 0. 



Then, equation ( [5.1| ) has the solution /(r) = k + cr"^, while ( |3.9D leaves N unconstrained, which is 
one of the A = oo solutions of 0. 

5.1.1 Associated Thermodynamics 



We abide to the assumptions of section to study the thermodynamics of the solution above. The 
relevant solution is the one with e = 1 with appropriate asymptotic behavior. The black hole horizon 
r+ is defined as the largest root of the equation /(r) =0, and therefore satisfies 



cr \ 



(5.13) 



The Hawking temperature is computed by the Gibbons- Hawking relation 



(5.14) 
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as 



27r 



4 + ^ 



1 + 



3?? / fc 



+ — 



(5.15) 



We may rewrite the mass as a function of the position of the horizon as 
M(r+) 





^3 




V 


3eG 





(5.16) 



We may now use the first law of thermodynamics to calculate the entropy by integrating dS = 
T~^dM, or S — So = fT~^j^dr^ . So is an additive constant to the entropy that we cannot fix from 
first principles. 

It is interesting to note that the entropies we will find below, whenever they accept the limit 
r+ — > (i.e. Y > 7^ 0), in this limit the entropies are finite. This is in contrast to various 

cases studied in the recent literature where the entropy diverges in this limit. For example, this is 
the behavior found in [|T^], that studied the entropy in the detailed balance and deformed detailed 
balance cases. We will therefore choose the additive constant below so that S{r+—>- 0) = 0. 

For k = (where we must impose rja > —2,^^ for the existence of the solution), integration gives 
directly 



G 



(5.17) 



For k = ±1 we have to distinguish the following cases 



l + ^>0 



G S{r^ 



Sn / k a 



3rjk 



:ln 



1 + ^- 



1 _LiZ2.. 



(5.18) 



1 + ^ < (where we must have ^>0) 



3ri\k\ 




1 + ^1 — + 



a 



(5.19) 



arctan ■ 



1 

J- + 252 1 



1 + ^ = (where we must have y >0) 

G S{r^ 



I Srjk 



6? 



2r . . 



(5.20) 
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5.2 3C + 87/ = 0, 



Having dealt with the case C + St] = 0, we may henceforth assume that C+St] ^ 0, and we will define 
the function gir^ as 

i .2 



fir) = k + 



4(C + 3r/) 



r + g{r). 



Equation ( p.4|) becomes 

The coefficients of equations ( |3.10| ) and ( |3.11D simplify to 



C+3r] 



Ai,=r[{3C + 8v)rg' + 2{C+A7])g] , Bn = {3C + 8r]) {r'g" -2g) . 
We now consider the next special case: 3( + 8ri = 0, Q-rj ^ 0. 

In this case, the non-linear term in ( p. 22 ) vanishes and the solution is easily found to be 



(5.21) 



(5.22) 



(5.23) 



g = r 



■ 2 / 6^- 



(5.24) 



where c is an integration constant. From ( p.23| ), we obtain Bn = 0, and therefore iV(r) = 1. A special 
case of this solution was found recently in |T^, while more special cases for this class of solutions 
were found in |jl4|. 

In order for the large distance limit to exist, we must have 



6^' 

T 



>o. 



(5.25) 



Expanding for r —>■ oo, using ( p.21| ), we obtain two asymptotic branches characterized by the sign 
e = ±1 



/(r) ~ + |^ei 
from which we obtain the mass as 



GM = - 



4 



2(C'^+6g2) 
3C2 



(5.26) 



(5.27) 



The solution with e = — 1 (c > 0) is defined for any r, while the one with e = 1 (c < 0) is bounded 
from below. In order for the solution to be asymptotically AdS when e = —1, we must have C < 0? 
and cr > 0. In the opposite case e = 1, we have C < 0, or ^ > and a > 0. 

When (7 = 0, one of the two solutions is asymptotically flat (which one depends on the sign of 
C). We obtain 

(c 



(5.28) 
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and matches the one found in HTll. 



For detailed-balance with A = 1, we obtain 



KdB^^^DB = ^ , OrflB + ^ = 0, (5.29) 



therefore the solution ( ^.21|) reduces to 

f{r) = k- Awr^ + Vcr , (5.30) 

which is one of the A = 1 solutions of [^]. 

A slight deviation from detailed-balance was defined in by 

C = Cdb(1-£') , v = VDB{l-e'), (5.31) 

where \e\ <^ 1, while keep their detailed-balance values. In this case, for A = 1, we still have 
3( + 8ri = 0, and indeed 

therefore, this slight deviation is a special case of the solution above. 



crnB + ^-^]= ^7^^>0, (5-32) 



5.2.1 Associated Thermodynamics 



The thermodynamics of the solution above has been studied in detail in several works [16|. We will 
therefore move forward to the generic case. 

6. The generic case 

Having disposed with the previous two special cases we may now consider the generic case in which 

(C+3r/)(3C + 8r/)^0. (6.1) 
We define the following combinations of the relevant parameters, 

(3C + 8r/)2 ' 3( + 8r7V2(C + 3r/) / ' 3C + 8r/ ' 4 ^ ^ ^ ' 



We define the function g{r) as before 



The differential equation (|5.22| ) now becomes 



r^g'^ + {C- A)rgg' + -(8 - C)g^ - Br'^ = 0. (6.4) 

In appendix y we describe in detail the derivation of the most general solution of ( |6.4|) as well 
as the one for the function A^. To describe the solutions we must distinguish two main cases and a 
few subcases. 
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6.1 A > 



In this case, it is -B > and < C < 6. Note that C = corresponds to the case examined already 
in section |5]1 . 

The general solution to equation (|6.4|) in our case is given in implicit form as 



7" \ 3 
To 



A 

1 — 

B 



C g 



2VB 



exp 



2eVA 
C 



arcsm \ 



^1 
B r2 



(6.5) 



where e = ±1 is a sign and > is the integration constant. The modified lapse function is 



N{r) 



C /r,x3 



3C 



1 Ag^ 



C 



PIA- 
arcsin-^' 3^ 



B r 



7^ 



(6.6) 

where we have chosen the multiplicative integration constant so that the function asymptotes to one 
at infinity distance. ( |6.5| ) in particular implies that \g{r)\ < \J^r'^- 

The plot of {r /toY against y = y for the e = 1 branch of ( |6.5|) is shown in figure |l], where we 



fixed C = 2^ A. Note that there are two branches. On the left hand side, r starts at rmin = e 
(where the curvature singularity lies) when y = —1, and increases till it becomes infinite at y = 
On the right hand side V > asymptotic infinity is at ?/ = and the radius decreases until a 
finite value r^m = c^Tq at ?/ = 1. For spherical or toroidal symmetry, k = 0, 1, and ( + Sr] > 0, the 
left-hand branch has a conventional horizon at which /(r) vanishes, for a range of couplings although 
the right-hand one does not. When + Sr] < it is the right-hand branch that has a conventional 
horizon for a range of couplings, while the left-hand branch has /(r) < 0. 



150 
125 
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. 5 



Figure 1: ^ (vertical axis) is plotted against U = (horizontal axis), for e = 1 and C - 

Changing the last ratio stretches the horizontal axis, but the topology of the graph remains intact. 



2VA. 
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Figure 3: On the left: the modified lapse function iV as a function of y for the left branch of the solution 
in figure 1. It asymptotes to one y = (asymptotic infinity) and diverges from positive values at the 

singularity y = —1. On the right: the modified lapse function as a function of y for the right branch of 
the solution in figure 1. It asymptotes to one y = (asymptotic infinity) and diverges from positive 
values at the singularity y = 1- Note that unlike the previous case it is no longer monotonic. 



The e = — 1 branch of ( |b.5| ) is shown in figure ^ for the same values of the parameters. It is just 
the mirror image of the e = 1 branch. For both signs of e, although the solutions start at a minimum 
Tmin, the curvature is singular there (y = ±1). 

In figure |^ we plot iV as a function of y for the solutions e = 1 above. At the curvature singularity 
y = ±1 it is iV ^ oo. 

For large distances r ^ (which makes sense for e{sgng) > 0), the leading behavior of (|6.5| ) for 
the solutions with positive mass is 



4(C + 3r7) V " 3 



4(T 



/(r) ~ A; + -—^ ^+ eWe-^(C + 3r?) r' + 0{r-^) , (6.7) 



2GM 



r 
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with 



Moreover we have 



2GM 



Brl (2^fA 
- exp 



'2A\ 



g = esgn{C + 3T]). 



N{r) = l + 0(r-^). 



(6.8) 
(6.9) 



This imphes that M = M. 
6.2 v4 < 

In this case, B can have either sign and C > 6 or C < 0. 

The general solution for g{r) in this case is given implicitly as 



^ .9^ Cg 



B-A^ + eV^^ 



9 



(6.10) 



where e = ±1 is a sign and r,, > is the integration constant. The modified lapse function is 



N{r) = 



\BC\ 
6 



\BC\ 



1 + 



c 



I) ^ 



\BC\ 



B-A^ + eV^^-i 



B~A 



B-A^-i-e 



cg_ 

2 r2 



B-A^ 



2 



(6.11) 



where we have chosen the multiplicative integration constant so that the function asymptotes to one 
at infinite distance. 

Plots of the solutions above for different values and signs of B and C are portrayed in figures ^ 
and 1^. In this case, r ranges from to infinity. There is also here a curvature singularity at r = 0. 




-10 -7.5 




(horizontal axis) for e = 1, ^ < 0, 

B > and ^ > 1. There are obviously two branches. The "cusp" corresponds to r = oo. The r = points 



Figure 4: On the left: r^/r^ (vertical axis) as a function of y 



(singularities) occur at y = ±oo. On the right, the same plot but with < < 1. Here, asymptotic 
infinity is at y = — oo, while the singularity is at ?/ = +00. 



For large distances r ^ r^, which for C > 6 makes sense for i? > 0, e{sgng) > 0, while for C < 0, 
B < 0, e{sgng) < 0, the leading behavior is achieved with the vanishing of the first factor in (|6.10|) 
and we obtain for the solutions with positive mass 
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Figure 5: On the left: r^/r^ (vertical axis) as a function of y = y (horizontal axis) for e = 1, A < 0, 

B > {) and > > — 1. Here, asymptotic infinity r = cxo is at y = cxd, while the singularity r = is at 
y = —CO. On the right, the same plot but with — ^ < —1. The "cusp" corresponds to asymptotic infinity 
r = CO. The r = points (singularities) occur at y = ±00. 



/(r) ~ + 



2GM = 



4(C + 3r7) 
3^ IBC 



4cr 

e+eVe-y(C + 3r7) 



2GM 



6 



1 + 



2J\A\ 



C 



VIA| 



e = es5(n(C + 3r7) 



iV(r) ~ 1 + 0{r-^) . 



(6.12) 



(6.13) 
(6.14) 



Note that to obtain the result ( |6.14|) one has to consider in the expansion of ( |6.12| ) also the exact 
term. 

This implies that also here M = M. 

For C < 0, there is another way to make something vanish on the left-hand side of ( |6.10|) and 
this is to have g/r"^ ^ 00. This gives a different asymptotic behavior at large r 



/(r) ~ A; + 



1 



4(C + 3r7) 



^^2 



(2a/|A|) r"^ 



|C| 



\C\ \ \C\-2^ 



3|C| 



+ 



N{r) ~ ( 



^ \C\+2^\A\ 



Note than in this case the exponent 



3|C| 



(6.15) 

is negative and spans all real negative values as — 00 < 

\C\-2^\A\ 

C < 0. This is a case where the long distance behavior in the /(r) function is stronger than the 
generic (9(r^) AdS behavior. 



6.3 5 = 

The value i? = is a special one for equation (|6.4|) with corresponding solution 



g[.r) = Cir 



2e^/\A\ +2 



C 



(6.16) 
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where e = ±1 and ci, No are integration constants. Note that i? = is allowed for the solution (|6.10|) 
but not for the solution ( |6.5|) . 

Detailed-balance with arbitrary A has 

and therefore corresponds to this last case. Substituting these values for A, C in ( |6.16| ) we obtain 

^(r) ~r(^^+"^^)/(^-^) , iv(r) ~ r-(^+^^+2£^^)/(^-^\ (6.18) 

which is the first solution found in with generic A > 1/3. 

When we choose further for the integration constant Ci = in (|6.16|) we find g{r) = 0, and then, 
for detailed-balance with arbitrary A, the lapse N remains undetermined, which is the second solution 
found in with generic A. 

We end by noting that in the generic case with the exception of i? = 0, there are always solutions 
with the standard asymptotic behavior exemplified in (|6.7|), ( |6.9| ) (or the same ( |6.12|) , ( |6.14| )). The 



solutions are generically asymptotically AdS/dS, with a standard 1/r tale that defines the mass of the 
solution. In all such cases the solutions come in a pair with different asymptotic effective cosmological 
constants 

3 



A 



4(C + 3r/) 



4cr 

^+hl^'-^iC + 5v) , e = ±l. (6.19) 



3 

Not only the bare cosmological constant is responsible for Ag// but also the curvature square couplings 
(, 1]. This is special example of what is known in cosmology as self-acceleration: the generation of 
acceleration (background curvature) by gravity alone. 

The effective cosmological constant can become zero and the solution asymptotically flat when 
(7 = 0. In this case, the e = — 1 solution is asymptotically flat, while the e = 1 solution has 

A.f = (6.20) 

2(C + 3r/) ^ ^ 

7. The horizons in the generic case 

We define again the exterior horizon position as the largest root of equation /(r) = 0. This implies 
that 

g{r+) k ^ 



r\ r\ 4(C + 3?7) 

where we defined the variable X that will be useful later on. 

7.1 A > 

From (^) we obtain an equation for r_|_ 



-X, (7.1) 



A C /r.N3 

I- X^ + e^X 
B 2VB 



(7.2) 
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The existence of the square root imphes an upper bound on r_|_ 



k 



< 



rl 4(C + 3r7) 

For a flat horizon A; = 0, the inequahty constraints the couphngs 

(1+3^) [2e-a{x+m>^- 

7.2 A < 

From ( |6.10 ) we obtain an equation for r+ 



(7.3) 



(7.4) 



^B-AX^ + e-X 
2 



WB-AX^ - ev^X|^ 



(7.5) 



8. Thermodynamics in the generic case 

We now investigate the relevant thermodynamics in the generic case following the remarks and as- 



sumptions of section |3^ . 
8.1 A > 

As showed before, there are solutions with the standard 1/r asymptotic behavior 
simplicity we assume that ^ + 3^] > 0, so that e = e. 

The Hawking temperature is computed from ([5.4[) as 



I). For 



I C-4 k 
8(C+3r)) ^ 2 71 



1-^X2 



[arcsiriY^ + e arcsin (y^^ X) ] 



(8.1) 



while the mass from 



IS 



2v^ V A 



[arcsiriY/p" +e arcsin(Y^X)] 



We calculate the entropy by assuming the first law dS = T ^dM, from which 



S = S,+ I T-'^dr+. 

dr^ 



(8.3) 



For k = 0, integration gives 



S 



6 A e 



C y 516(C + 3r7)2 G ' 
where we have set 5*0 = so that S'(r+ ^ 0) = 0. 



•4) 
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For k ^ (where the inequahty (|7.3|) is relevant) 



S{r+)-So 



■Kk QA 



G \ BC 



kr , 



B 



/ A " 
— + arcsinf \ — X 
\\ B , 



+ 



A 



■ arctan(h) 



4(C+3r?) B_ k 



A ^ 



where 



a = e- 



4(C + 3r/) /S 



while arctan, arctanh are used for |a| < 1, |a| > 1 respectively. 
8.2 A < 



(8.6) 



Also here, there are solutions with the standard 1/r asymptotic behavior (|6.12| ),( p.l4| ). We will study 
the case with C > 6 and for simplicity we assume that C + > 0, so that e = e. 
The Hawking temperature in this case is computed to be 



T(r+) 



4:71 



|_ C-4 k 



8(C+3r)) 



2 ri 



1-^X2 



+ eVB 



X2-f -eX 



'X2-f + eX. 



•7) 



while the mass is 



M(r^ 



No M 

G Vec 



2vPT V A 



X2-f -eX' 



'X2-f + eX. 



Using the first law, we calculate the entropy by integration as before. For k = 0, this integration 
gives again ( ^^) . 

For k 0, we obtain 



S-S. = — 



nk Q\A\ 



G \ BG 



krUIX^-j 



elog 



'X2-f + e(a-v/TTi?)X-J|^ 



'X2-f + e(a+v/T+^)X-J||| 



where 



a = e- 



_4(C + 3r7) B 



\A\ 



We note that the expression ( p.8[ ) has a logarithmic divergence when r_| 
to the first factor inside the logarithm. 



for both branches e, due 
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9. Discussion and further directions 



We have found the most general spherically symmetric solution to Hofava-Lifshitz type of gravity 
theories with zero shift. The action used was general except at the level cubic in curvatures where 
the Cotton^ combination was used. We do not expect new qualitatively different effects if we allow 
also the most general (curvature)^ couplings. 

• It is a generic feature of the solutions found that they have regular large distance asymptotics 
that are asymptotically AdS/dS or flat. Moreover, generically the next correction is compatible with 
a standard Newton's law. 

There are exceptions to this results that are special. In particular, the detailed balance action 
first written down by Horava is one of these notable exceptions. As shown already first in ||^ it 
does not reproduce the correct Newton law at large distances. This is part of a class of special cases 
analyzed here, occurring when 



although in the second case it is not the generic solution. 

• In one of the two general categories of parameters {A < 0) where the horizon distance is 
permitted to shrink to zero, the entropy has a logarithmic divergence in this limit. In the other 
general category {A > 0) the horizon distance is bounded from below. 

• For special values of the curvature- squared couplings ( + 3i] = 0, when the horizon position is 
allowed to go to zero the entropy is regular in this limit. 

• For the same values of the curvature-squared couplings C + = 0, (-rj ^ 0, there exists a 
solution with logarithmically corrected large distance asymptotics. Its "mass parameter" defined in 
the naive way, depends logarithmically on radial distance and becoming smaller at larger distances. 

• When the cosmological constant term in the action is absent, the effective cosmological constant 
of {^AL)dS curvature is also zero (for some branches). 

• We have studied geodesies of particles with finite and infinite light speeds in spherically sym- 
metric backgrounds with traditional horizons. We parameterize the dispersion relations of particles as 
Po ~ (P^)") with n > 1. We find that when n> \ the traditional behavior of the horizon disappears 
suggesting that for such particles that black hole is effectively "naked" . 

There are several interesting questions remaining open. 

• The meaning of the notion of a black-hole and that of horizon in a theory with particles that have 
modified dispersion relations must be investigated. In particular, apparent puzzles/instabilities 
that seem to exist in such cases must be resolved or such theories will be physically discredited. 
A key element that may turn out to be important is a strong downgrading of energy in cases 
of modified dispersion relations that together with low energy Lorentz invariance may be a way 
to resolve puzzles. 

• Finer properties of the black hole found should be studied as well as extremal limit in search 
for stronger constraints on the physical acceptability of such solutions 




or C 



16(C + 3r7) 
3C + 8r7 



< 0. 



(9.1) 



- 24 - 



Spherically symmetric solutions in the projectable theory need to be studied as no non-trivial 
one is known so far. 



The study of the non-linear structure of spherically symmetric solutions in the modified non- 



projectable theory of must be studied in order to evaluate the stability claims of such 
solutions at the non-linear level. 



We plan to return to these question in the near future. 
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APPENDIX 

A. The classical equations of motion 



We now add the action of matter 



Sm = J d^xdty^N C^,^,ter{N,Ni,g,j) (A.l) 



to the gravitational action ( |2.11| ) and we will vary with respect to the gravitational fields to obtain 
the equations of motion. 

The equation obtained by varying N is 

-a {KijK'^ - \K^) + pCijC'^ + -i£'^^RiiVjR\ + CRijR'^ +7]R'^ + ^R + a = Jn, (A.2) 

with 

Jn = ^-^matter ~ N — — . (^-3) 

OiV 

The equation obtained by varying Ni is 

2a(VjK^' - XV'K) + at^^^ = o. (A.4) 

oNi 

Finally, the equation of motion obtained by varying gij is more voluminous 



1 



(^S Qmi m)',kin9j£ Qmi)-!^^ 9jf- Qmi) -jkRji 



— {S"^'^^QmiRk)-ngjl + {£"^^^Q'^Rki ) ■,ngje + - {^'^''^R\keQm ) " QklC'^Qij + 

£""''QmiRjE;k] +a[N{2r]R+0]9ij + N{2r]R+ORij + 2N{CRikR--PCikC!;) 
-[Ni2r]R+Oh + a[Nmj + lQ,)]-2[Nmk + lQk)]/ + [NiC^^^ 

-^{(3CkiC'''+^Rk,C''' + CRMR'''+vR^+^R+cy)g,j+2QN{K,kK';-\KK,^) 

''^iKkeK''-XK')g,, + ^g,kg,E^\^{K''-XKg'')] + a[{K,k-XKg,k)Njf 



2 ' 

+a[{K,k-\Kg,k)N,\^-a[{K,,-\Kg,^)Nkt + {i ^ j) = -2N^-^^, (A.5) 

where S^^^ was defined below equation ( |2.5| ), ; stands for covariant differentiation with respect to the 
metric gij, and 

Q,j=N{^Rij+2pCij). (A.6) 
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B. Derivation of solution when C+3r] = 0, (-t] ^ 0. 
In this case, equation ( |3.4| ) becomes 



r""/'^ + A\-r^-(f-k)]rf' + A\^r'' + (f-k)](f-k) - —r^ = 0. 

IT] 1 iT] J 1] 



We define a new function y and a new coordinate R 



Equation ([B.l|) becomes 



fir) = k -\ H — r'^y{r) , = Inr . 

6C V 



?/2 + 4y+12y = 0, 



where a dot is a derivative with respect to R. 

The general solution of equation ( |B.3| ) must have y(r) < 1/3, and has two branches with e 
It is given by the implicit expression 



r 



T \ 3 



where is integration constant (positive or negative). 

We proceed with the integration of equation (|3l9|) and define 



6e 2e V r 



y 2e(yr^-e) 



where in the last equality we used equation (|B.4D . 

Moreover, using (|B.3[), we obtain y = — 4F(1 + 3Y), and therefore 



Equation ( p.9|) now becomes 



Equation ( [B.7| ) can be integrated to 



Y 



dlnN 
dY 



N-- 



1 + 3F 
l + QY' 

-123+*' 



1 + 6F 



No 



e^^\l + QYy 

with A^o integration constant. In terms of the y variable it becomes, 

pf(v'T^-l) 

N{r) 



yi^ 

where we fixed the integration constant = e~^~^. 



Equation (p.9| ), together with ( p.4|) provide the general solution when ( + 3ri = 0. This solution 
can be given alternatively in terms of the variable Y as 

f(r) = k + —r^ -—r^Y{l + 3Y) (B.IO) 

-^y^^_(!jL)' N(r) = ^ 

2\rJ ' ^' e^^\l + QY\ 



l + 3F)e-^^ = - - , Nir) = —-^. (B.ll) 



Above, the branch with e = 1 corresponds to F < —1/6, while e = —1 corresponds to F > —1/6. 
We will now study the large distance, r —>■ oo behavior of the solutions above. 

• For large distances r » \ro\, the e = 1 branch of (|B.4| ) can be approximated as ?/ ~ r~^, and 
therefore 

/(r) ^ + ^r^ - ^ + 0(r-=^) , GM = (B.12) 

For an asymptotically AdS or flat solution we must therefore impose a > 0. 
For the modified lapse and the lapse we obtain 

iV(r) ~ 1 + C(r-^) , Ar2(r) ~ A; + — r2-^^ + 0(r~3). (B.13) 

6^ r 

• In the case with e = — 1, the large distance limit r — > cxd is achieved when y — > — oo. We find 
to leading order that 



^r^^log(^) -loglogf;^) +■■■ , iV(r)^ (-)" + ••• (B-14) 

W VoJ |log^(r/ro)| V r / 

Therefore, such terms logarithmically dominate the cosmological constant in the large distance 
limit. 



To obtain the entropy expressions ( |5.18| )-( |5?^ ) of the above solution for /c 7^ 0, it is helpful to 
convert first to Y_^ 

S-S^= fr^T-'.^^dY^, (B.15) 



where we can find 

vre^r^ e6^+(l + 6y: 
~6r2+5y^ + l- 

1 k ,k, a,(l + 6F+)(6y| + 5F+ + l-J^ 

so finally. 



kr]?!- fk a\ f (1 + 6F+)^ /t^-,on 
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C. The general solution 

In this appendix we describe first the general solution of equation (|5.22| ) that we reproduce here for 
convenience. 

(3C + 8r7)rV' + A{C + Ai^)rgg' + A{C + 2i^)g^ + ^ {Aa-^^r^ = 0. (C.l) 
We change variables as 

2{C+4»;) r. 2{C+477) 

h{R) = e^<+^^g{R) =r^^c+Wg , i? = lnr, (C.2) 
with dot being the derivative with respect to R. We also introduce 

A , «<K+3.rt ^ 1 ( 3^1^ X ^ 16(C+3„) ^ac(6-C), (C,3) 



so that equation ( |U.1| ) is 



rV + iC~ 4)rgg' + ^(8 - C)/ - Br' = 0. (C.4) 



Under ( |C.2|) equation (|C.1| ) becomes 

/i' + Ah^ - Be^^ = 0. (C.5) 

From (|C.5|) we observe that Be^^ — Ah'^ = Br^ — AK^ > 0. If W{R) satisfies the first order 
differential equation 



2 

then, /i(_R) is given as 



CR 

This is the general solution for h{R). This is shown by defining 

mn) ^ (C.8) 

and showing that equation ( p.5|) implies that W{R) satisfies equation ( |C.6|) , and h{R) is given by 

Since (|C.7D gives h'^, only |/;,| will be determined from the solution, and therefore, two branches 
will always arise for each solution \h\. From ( |C7| ), we must have B /{A + W""^) > 0. 
Equation ( |C.6|) leads to a quadrature 

dW r 
c .TTTT— ^777^ = ^-^0 = In-, (C.9) 



1-^W){1+AW^) To 
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where > 0, Ro are integration constants. The integrant can be broken as 



1 



1 



2ACW 



+ 



AA \ 



{1-^W){1+AW^) AA+C^\1-^W l+AW^ l+AW. 
The integration of the first two terms in ( p.lO| ) gives 

\l + AW'^\ 



(C.IO) 



6 



In 



(C.ll) 



The integration of the third term depends on the sign of A. There are two cases, one for A > 
and another for A < 0. 

Before we proceed to consider the two cases, we should also consider the second equation, namely 
that in terms of the current variables, its coefficients translate as 



An = (3C + 8r/) 



rg_ 
W 



Therefore, equation (|3.9|) is rewritten as 



Bn=2(3C+8„)^(l-|iy" 



2A l-^W 



(C.12) 



(C.13) 



dW C 1 + AW^ 

We now consider the two cases separately. 

C.l A > 

Equation ( |C3| ) implies that < C < 6 and a solution exists only for i? > 0. The solution for W is 

'ro\6 1 + AW^ 



In 



rJ (l-^W) 



+ arctan (/A W) = 0, 



(C.14) 



and for h{r) we have the implicit equation 
In 



^ A h'^ _^ C h 



2eVA 

LVTo/ V B r'^ ^2\fBr'^l'^ -I C 
with h taking both positive and negative values and e = ±1. 

2(C+4>?) _ C 
3C+8»? 2 

^2t^fA 



arcsm 



'A h 



B 



(C.15) 



Using ^yts!!^ ~ rewrite the solution (|C.15| ) as an implicit equation for g{r) as 



^ \ 3 



B 2\^Br 
Integration of equation ( |U.13| ) gives 

N{r 



exp 



A g 



(C.16) 



In- 



arctan(vy4iy) 



(C.17) 
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and using (|C.14|) we obtain 
and in terms of g{r) 



iV(r) = No 



\i-qw\ 



(C.18) 



^ Br* 

Finally, we choose for convenience the constant of integration as 



1 Ag2 ^ C g 



Br* ^2VBi 



No- 



e c 



arcsin(^ B 



1 A 52 

-L n „d 



(C.19) 



so that the final solution becomes 



exp 



^{sgng)—^ arcsmW — ^ 



(C.20) 



6 V r . 



e C arcsmy ^ 



1 ^A31 
^ B r* 



1 Ag-^ _ ^ C g_ 



(7 g-^^[arcsin(^^)-(ssng)arcsin,/|g"] 

6" ^ 



1 _ £121 



(C.21) 

Equations ( |(J.16D , ( |C21| ) form the general solution of the system for A> Q. 
For large distance r » (which makes sense for e{sgng) > 0), the leading behavior of ( p.l6|) 
for the solutions with positive mass is 



9 



2B s/Brl (l^fA 



-2 ~ 3C 3r3 



exp 



2A\ 




We therefore obtain 



4(C + 3r/) L 



^+eVe-y(C + 3r7) 



with 



2GM 



exp 



2VA 




. 2GM 4, 

r 



, e = es5(n(C + 3?7). 



Moreover we have 



iV(r) = 1 + C»(r-^). 

To obtain the entropy expression (|8.5|) for A; 7^ it is helpful to convert to W^, 

1 dM 



S - So= / r+T 



(C.22) 

(C.23) 

(C.24) 
(C.25) 

(C.26) 



where we can find 



r,T-^ = ATT 



L 2(C + 377) 



'5(1- 



C-4 



^ e f arctan( 1^+ ) -arcsin^2A/3C ] 



(C.27) 
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1 dM 



B k 



r+dW+ AGl+AWl 
so finally, 



2^ r 



2(C + 3r/) 



i^2A/3C-arctan(v^Vy+)] 



1 2 



(C.28) 



4(C + 3r/) 



G y 3C J 1+AWl 

This integration can be done explicitly in terms of defining tan^ 
terms of r+. 

C.2 A < 



(C.29) 
A W+ and replacing in 



\1+AW^ 



and in terms of h{r) 



B-A— - e- 



2 r<^/2 



/l2 



IS 






6 


(a 




- h 


c 


\c/2 





T \ 3 

r 



where e = ±1, ro > is the integration constant, and = \B\2~ 
Finally, the solution ( P.31|) is written in terms of g{r) as 



B-A^ - e-^ 



2 



r"* 2 

where g is allowed to take both positive and negative values 
Integration of equation ( |C.13|) gives 



7^ \ 3 



^/\1 + AW^ 



1 + 



and using ( |C.30|) we obtain 



N{r) 

which can be rewritten in terms of g{r) as 
N{r) 



3 1 + AW^ 



c 



No 

Choosing No 

Mr) = f^^' 
V r 



T \ 3 

r 



\B\ 






^ r-i 2 r2 



1 + 



2^1 



c , we obtam 



\BC\ 



\BC\ 



\BC\ 



B-A^ + eV^^i 



B-A^ 



2 

7^ 



1 + ^1)^ 



(C.30) 



(C.31) 



(C.32) 



(C.33) 



(C.34) 



(C.35) 



(C.36) 
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Equations ( C.32 ), ( C.3(: ) form the general solution of the system for A < 0. 



For large distance r»ro (for C > 6, i? > or C < 0, i? < 0) the leading behavior of ( |C.32| ) for 
the solutions with positive mass is 



g[r) ~ e 



3C 

We therefore obtain 
/(r) ~ A; + 



~ — + 0(r~^) , 2^714^ = — 

r 3 



1 + 



2^/\A\ 



c 



e+eVe-y(c+3r7) 



4(C + 3r?) 
Moreover we have 

iV(r) = l + 0{r-^). 

To obtain the entropy expression (|8.8|) for A; 7^ it is helpful to convert to W+, 

1 dM 



-dW+, 



where we can find 



4tt 



+ eVB 1 



C-4 



1 rfM No SB 



r+dW+ 20"^ 2C l + AWl 
so finally, 



^ ^ 37rA; /25 
b - bo = e—-\ — 



(C.37) 



+ C(r-^) , e = esc/n(C + 3r/). (C.38) 



(C.39) 



(C.40) 



2 ' 



(C.41) 



(C.42) 



(C.43) 



G \ 3C J 1+AWl 

This integration can be done explicitly in terms of W+ defining u = tan(6'/2), sin^^ = \^—A W+ and 
replacing in terms of r+. 



D. Dispersive geodesies 

Particles with non-standard dispersion relations, and in particular non-Lorentz invariant ones, do not 
follow the usual geodesies in the gravitational field. This has been discussed in several very recent 
works ||26|, |2^. We will follow a slightly different approach here to study dispersive geodesies. 



We will eventually investigate the fate of the horizon for such particles and find that it seems to 
disappear, as intuition suggests. 

To proceed, we write the metric in the ADM form 

Goo = -N^ + N,g'm, , Go^ = N, , G,, = g^, , det[G] = detf^^jiV^ (D.l) 
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as well as the inverse metric 



^00 _ ^ 



= g'- 



(D.2) 



where indices of Ni,gij are raised with the spatial metric gij. 

We will also consider a scalar mode which in the standard relativistic case has a coupling to 
gravity in ADM notation given by 



S = / d'xdt^N 



(D.3) 



For the non-relativistic case we write a general non-relativistic quadratic action in the Lifshitz spirit 
as 



d^xdt^N 



1 

'iV2 



, U = g'=VNj 



with equations of motion 



H $ 



-A 



\/9 

N 



($ - N'di<!>) 



-FP$ = 



(D.4) 



(D.5) 



where in the spatial part we neglected derivatives acting on gij,N as they will be irrelevant in the 
derivation of the geodesic (geometric optics) equations. 

To pass to point-like trajectories, we replace idt po, idi Pi, and neglect metric derivatives 
to obtain the equivalent Hamiltonian ( "zero energy" ) constraint 



H 



{po - N'p 

iV2 



+ F[C]=0 , ( = g''p^p, 



This is implemented with a Lagrange multiplier e, to obtain the world-line action as 

»i 



dr 







Pot + pix' + -H 



(D.6) 



(D.7) 



where r is the affine time of the path and dot stands for dr- e is a one-dimensional einbein. The 
action ( p.7| ) is invariant under general reparametrization of the affine time 



r = f{r) 



fir) 



(D.8) 



Using such reparametrizations the einbein can be set to a constant T. 

The classical path equations come from varying e, t, x*, and po, Pi. The variation of e by definition 
gives ( p.6| ) while the momentum variations give 



. e6H . e , 







eSH eN' 



. 

Po - N'Pi = —t (D.9) 
+ NH = -eF'(C) g'^Pj (D.IO) 
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Finally the variations over the coordinates give 

edH 



Po 



2 dt 



Pi 



edH 
2dx^ 



Equation (p.lO|) gives 



C(i^'(C)) 



2 g^j{x' + NH){x^ + NH) _ ^ 



(D.ll) 



(D.12) 



This is a non-linear equation in general, and we can pick a solution C(^/e^)- Then we can express 

g,,{x^ + NH) 



Pi 

We also obtain by substitution 

iV2 



N"^ . NJx' + NH) 



(D.13) 



Pot + p,x' = —e - ^ ^ = —t^ - eCF'iO 



(D.14) 



where in the last equation we used ( p.l2|) . We also have 



H = -—t' + F{C) 



(D.15) 



and finally the on-shell action is 



wl 



dr 



P + e(F(C)-2CF'(C)) 



(D.16) 



where C, should be thought as a function of the metric, i;% and the einbein as given from ( p.l2|) . 
Indeed varying e,x\t using the constraint (p?T|) gives that appropriate geodesic equations and the 
zero energy constraint. This action reduces to the standard relativistic action for F{() = C + 

We now consider a particle with a non-standard dispersion relation F{() = (corresponds to 
a dispersion relation {p^Y — (p^)" = 0) a radial geodesic and a spherical symmetric BH metric to 
obtain 



S=- f dr 
2 Jo 



NH 



212 



+ (1 — 2n)e 



2n^ 



2N 



2n^ 

where E' is a constant of integration (energy). 
From (p.2|) we also obtain 



. 2 f2N 



2n^ 



(D.17) 
(D.18) 



(D.19) 



For any n > 1 the geodesic is regular at the zeros of /, N . 



- 35 - 



If at r* / has an a-th order zero and A^^ a b-th order zero, then the geodesic is singular if 

a + ->2. (D.20) 
n 

We have always b > a. In the case of HL we expect generically that n — 3, therefore the condition is 
a + b/3 > 2. For a single zero of /, a = 1 and then 6 > 3 that is unhkely. Generically b — a — 1. 

For a double zero of /, then we always have a singularity. This suggests that only extremal BH 
have a horizon in support of the intuitive notion that if a particle has a dispersion relation with an 
infinite speed of propagation, then it feels no horizon in the background of BH solutions. 
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